NACA TN 2501 gggg 


NATIONAL ADVISORY COMMITTEE 
FOR AERONAUTICS 


TECHNICAL NOTE 2501 


EXPRESSIONS FOR MEASURING THE ACCURACY OF APPROXIMATE 

SOLUTIONS TO COMPRESSIBLE FLOW THROUGH CASCADES 

OF BLADES WITH EXAMPLES OF USE 

By John T. Sinnette, Jr., George R. Costello and Robert L. Cummings 

Lewis Flight Propulsion Laboratory 
Cleveland, Ohio 


NACA 


Washington 

October 1951 AFMDC 

TECHNICAL LIBRARY 
AFL 2811 


I 


NATIOML AWISOET CCMMiraEE FOR AEROMUTICS 


TECH LIBRARY KAFB, NM 




TECHNICAL NOTE 2501 


EXH?ESSIONS FOR MEASURING THE ACCURACY OF APPROXIMATE 
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SUMMARY 

Four necessary conditions for steady irrotational ccurpressihle flov 
throu^ a straight cascade of blades are derived from the irrotationality 
condition and the conservation of mass and momentum. Expressions are 
obtained which measure the degree to which an approximate compressible- 
flow solution departs from these conditions. The expressions may be 
used as a measure not only of how well a given flow solution approxi- 
mates the flow of the usual polytropic gas^ but also of how well it 
approximates the flow of an arbitrary barotropic fluid. 

As Illustrations of the use of the error expressions, they are 
applied to three basically different types of approximation to the flow 
of a polytropic gas through a cascade of typical compressor blading; 
namely, (l) the incompressible approximation, (2) the Prandtl-Glauert 
approximation, and (3) an approximation based on a linear relation 
between pressure and specific volume. The approximation based on the 
linear pres sure -volume relation gave mich the best agreement by satis- 
fying the Irrotationality and continuity conditions exactly and the • 
momentum conditions within about 1 percent. Although it is shown that 
in the Prandtl-Glauert approximation, the four necessary conditions for 
the compressible flow are satisfied to the extent that terms higher than 
the first order in the perturbation velocity are negligible, the errors 
resulting from this linear approximation were so large for a typical . 
compressor cascade that this approximation was no better than the inccm- 
pressible approximation (maximum error ejqjression, 19 percent) . When 
both parallel and normal components of the perturbation velocity were 
considered in computing the resultant campressible velocity, however, 
appreciably better accuracy (maximum error expression, 7 percent) was 
obtained. 
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INTROIXJCTION 

The c empress ible flow along the blade siirface must be accurately 
determined In designing blade rows for high loading and high subsonic 
Mach numbers in order to avoid excessive local velocities and to prevent 
flow separation. Because the solution for the two-dimensional potential 
flow of a compressible perfect fluid through a cascade of blades cannot 
be obtained analytically, some approximate method must be used. The 
methods of approximation consist in replacing the original partial dif- 
ferential equation by either a simpler partial differential equation or 
a difference equation and solving. !Bae solution of the simtpler equation 
may Involve additional approximations , such as modified boundary condi- 
tions. 

Some methods of approximation give the flow throughout the field 
and others permit the detennlnation of the flow on the blade surface 
directly. For all methods, however, it is desirable to obtain high 
accuracy on the blade surface where the conditions are the most critical 
because there the velocity reaches its maximum value and its distribu- 
tion along the blade surface determines the boundary-layer behavior. In 
general, the errors involved in the approximate methods are of two 
types: (l) errors arising from approximating the exact-differential 

equatlonj and (2) errors arising in the solution of the simplified equa- 
tion. Because the approximate solutions do not usually give in them- 
selves a convenient means of checking the over-all accurapy obtained, an 
investigation was conducted at the KACA Lewis laboratory to obtain cri- 
terlons for measirrlng and comparing the accuracy or consistency of the 
various approximate contpress ible -flow solutions. In a manner similar to 
that used by Weinig in reference 1 to obtain checks for the numerical 
accuracy of incompress ible -flow solutions, four Integral relations based 
on the necessary conditions of lirotationality and conservation of mass 
and momentum were obtained. Error expressions which give the deviation 
of an approximate compressible-flow solution from these necessary con- 
ditions are presented for meas'urlng the accuracy of campresslble-flow 
solutions. Because the expressions are given in a general form which 
permits using an arbitrary relation between pressure and density, they 
may be used to estimate how closely a given solution approximates the 
flow of any arbitrary barotroplc fluid. 

As an illustration of the application, the Integral expressions are 
used to evaluate the relative accuracy of three different methods of 
approximating the usual compressible flow of a polytropic gas (a perfect 
gas with constant specific heat) through a cascade of blades. The three 
methods of approximation compared are: (l) the incompressible potential 

flow approximation, (2) the Erandtl-Glauert approximation as adapted to 
cascades by Woolard (reference 2), and (3) an approximation based on the 
linear pressure-volume relation as applied to cascades in reference 3. 
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THBOET • 

The relations for conseorvation of mass and momentum, for steady two- 
dimensional coHrpressihle flow may he written 



(-pV^ dy + pVy dx) 


= 0 


( 1 ) 



p dx + pVyC-Vj^ ^ + "'^y 


= 0 


( 2 ) 


and 



dy + pV^(-V^ ^ + Vy dx)_ 


= 0 


(3) 


where the contour integrals are evaluated around any circuit 0 enclos- 
ing only fluid, that is, not enclosing any blades. (All symbols are 
defined in appendix A. ) The irrotatlonallty condition, which requires 
the circulation around the path C to he zero, may he written in the 
fonn 


(V^ dx + Vy dy) =0 (4) 

Equations (l) to (4) would, of course, he true for either direction of 
integration along C, hut this direction was taken as clockwise for the 
sake of definiteness. 

If the contour Integrals throu^ a cascade of blades are evaluated 
along the contour C Indicated by arrows along ahcdefgha in figure 1, 
the portions of the integrals on the free sections of one stagnation 
streamline ha and he cancel the portions on the free sections of the 
other stagnation streamline fg and de because the Integrands are 
equal at corresponding points and the integrations are in opposite 
directions. The only remaining portions of the contour integrals are 
those along ah and ef, which give the Integrals around the blade, 
and those along cd and gh, which can he readily evaluated in terms 
of the uniform upstream and downstream conditions. The contour integrals 
in equations (l) to (4) thus become 
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oV.. iv = foV..V 


( 5 ) 



p dx + pVy(-V^ dy + Vy. dx) 



dx + 





pV V 
y X 


<3y 



P dx + (pVyV^)^S - (pVyV^)2S 


( 6 ) 




+ PlS - PgS + (pv/)iS - (pV^2)^s 


(7) 



(Vx dx + Vy dy) 






7i/ 


( 8 ) 
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where B indicates the integration path around the hlade contour 
(counterclodcwise) and Vg is taken as the conrponent of the velocity in 
the direction of Integration or increasing s. 

For any potential flow in which the density is a known function of 
the pressure, the density and the pressure can he expressed as a func- 
tion of the velocity hy means of Bernoulli's equation 



(9) 


Consequently, if the approximate compressible-flow solution is given in 
terms of the velocity distribution, all the separate portions of the 
conto-ur Integrals can be evaluated and the magnitude of the deviation 
from the necessary conditions (l) to (4) determined. These deviations 
are expressed as follows: the percentage error in mass flow 


8m 


Wz - (pVi 


X 100 


( 10 ) 



and 



( 12 ) 


and the percentage error in circulation 
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Op = 



ds - (Vi^y 


'2,y 


)S 




)S 


X 100 


(13) 


EVALUATION CF TEE?EE TYPES OF COMPRESS IBEE -FLOW 
APPROXIMATION USING ERROR EXPRESSIONS 

Three different methods vere used for approximating the compres- 
sible flow through a cascade of typical compressor blades to serve as 
an illustration of practical application of the error expressions (lO) 
to (13) and their relative accuracy was estimated by the use of these 
four error expressions. 

IncQgpressible flow as approximation to compressible flow . - For 
many applications^ incompressible-flow solutions have proved to be 
satisfactory approximations to the desired compressible -flow solutions. 
Because of this observation and because of the relative simplicity of 
inconrpresslble flow, until recently practically all the extensive 
literature on cojiputed flow through cascade’ has been^evoted to incom- 
pressible flow (most of the literature prior to 1949 is siMmarized in 
reference 4). The first approximation method investigated is therefore 
an incompressible approximation. In this approximation, the dimension- 
less compressible velocity V^/ap is assumed to be egual everywhere to 

the velocity Vj_ for some known Inconpressible-flow solution. Inas- 
much as the velocity for any inccmipresslble flow may bre altered by a 
scale factor, a family of ccmpressible-Tlbw approximations can be 
obtained from any known Incompressible solution. The accuracy of the 
approximation decreases as the scale of the dimensionless compressible 
velocity V^j/ap, and hence the Mach number M, is increased. 

In order that the ex^ple illustrating the application of the ~ 
error expressions might be of practical Interest, the cascade and the 
scale of the velocities were chosen to be representative of axial-flow 
compressors. The incompressible flow was obtained by a numerical solu- 
tion using the theoretically exact Inverse method of reference 5 wherein 
the blade shape is coirputed from the prescribed velocity distribution. 
The velocity distribution chosen and the resulting cascade geometry for 
this example are shown in figures 2 and 3, respectively. In this 
example, as well as in those following in this report, the total blade 
arc length was taken as 2jt. The maximum dimensionless compressible 
velocity on the blade obtained by this approximation is 0.7558, which 
corresponds to a Mach number of 0.8031. 
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In order to oltain measures of the accuracy of this approximation 
to the usual compressihle flo'w of a polytropic gas, the error exj)res- 
sions (lO) to (13) were computed using the relations for the adiabatic 
flow 


and 



(14) 



(15) 


These relations may he readily obtained from the general Bernoulli equa- 
tion (9) using the adiabatic relation for a polytropic gas 



(16) 


Because the blade-profile velocities of the incompressible flow on 
which this approximation was based were computed at equally spaced 
points on the unit circle into which the blade profile is transformed, 
the integrals also were expressed in terms of the unit circle in order 
that Simpson’s rule might readily be applied. Thus 
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Although, the incompressible solution on which this approximation to 
conpressible flow was based was obtained by a theoretically exact method^ 
there still remain numerical inaccuracies due to rounded-off values, 
approximation methods for numerical integration, etc. It-was therefore 
desirable to have some measure of the accuracy of the Incompressible 
solution Itself. The accuracy of the incompressible solution was meas- 
ured by the same expressions (lO) to (13) as for the compressible 
approximation but ^riLth constant density and with the pressures obtained 
from Bernoulli's equation for incompressible flow 


P = Pt - ^ (18) 

The numerical values of the error expressions obtained for the 
Incorrpresslble solution and the incompressible approximation to a com- 
pressible flow with r = l-i are shown in the following table: 


Method 

Error 

(percent) 



^x 

6p 

Incompressible flow 

0.00 

-0.09 

- 0.06 

-0.07 

Incompressible approxi- 





matlon to compressible 





flow 

6.72 

2.05 

-14.14 

-0.07 


A comparison of the relative errors in this table indicates that the 
numerlceil inaccuracies in obtaining the incompressible solution are 
quite insignificant compared with the Inaccuracies Involved in this 
approximation to compressible flow. The large percentage error in the 
x-component of force Indicates that the inccmpresslble flow is a poor 
approximation to compressible flow through the cascade at this Mach 
number. 

Prandtl-Glauert approximation . - The incompressible solution was 
next used to obtain a compressible-flow solution based on the Erandtl- 
Glauert approximation as applied to cascades by Woolard (reference 2) . 
This approximation is based on the linearized equation for the 
pert\rrbatlon-veloclty potential and uses only linear terms in the 
perturbation-velocity components. Woolard considers the vector mean of 
the upstream and downstream incompressible velocities Vjj^ ^ as the 

reference velocity and deviation from this velocity as the perturbation 
velocity. If the transformation 
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is used where and X is the coordinate parallel to, 

and Y, the coordinate normal to the vector mean velocity, a linear 
w approximation to a compressible flow throu^ a cascade is obtained from 

w an inconpressible flow through a related cascade. The compressible cas- 

cade has the same chord and spacing in the X-direction Sj as the 
incompressible cascade, but the spacing in the Y-direction Sy is 
increased by the factor l/Q. The ratio of the compressible to the 
inconpresslble values of blade thickness, camber, and air turning angle 
will be equal (within the linear approximation) to 0/k. It is shown 
in appendix B that the four Integral expressions (l) to (4) are satis- 
^ fled within the linear terms in the perturbation-velocity components for 

the preceding Interpretation of the Erandtl-Glauert rule. 

' For the linear approximation to be reasonably accurate, however, it 

is necessary that the angle of attack, turning angle, blade camber, 
blade thickness, blade surface slope, etc. be sxifficlently small that 
the conpresslble and incompressible perturbation velocities be very 
small compared with the mean velocity Vm- In most practical cases 
these conditions will not be strictly satisfied everywhere and, there- 
fore, expressions (lO) to (13) serve as convenient measures of resulting 
mean errors in any particular example. As an illustration of this use, 
the preceding interpretation of the Prandtl-Glauert rule for the case 
of k =0 was applied to the IncompreBsible-flow solution discussed 
previously, and the error expressions were determined for this approxi- 
mation to compressible flow. In this case (k =f5), the blade section 
and the turning angle are the same (within the linear approximation) as 
for the incompressible flow and for the incompreBSlble approximation, 
but the blade spacing Sy is increased by the factor 1/n. Because of 
this change in Sy, the flow angles with respect to the normal to the 
cascade axis are also changed. The relation between the mean flow 
angles is given by 

■fcan tan 

The mean Mach nuniber was taken the same as for the Incampresslble 
approximation to compressible flow in order that the error expressions 
might be used to measure the relative accuracy of the two approximations 
under comparable conditions. The cascade geometiy for the compressible 
, flow is shown in figure 4. The blade spacing and orientation for the 

Incompressible cascade are also shown for ccmparison. 
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Inasmuch as the Prandtl-Glauert theory is based on the assumption 
that the pertiirbatlon velocity is sufficiently small that powers hl^er 
than the first in the perturbation-velocity components may be neglected 
throu^out, the use of other than linearized expressions in these ccm- 
ponents in computing the resultant velocity might seem unjustifiable. 
For most practical cases, however, quite different results are obtained 
if the linearized expression . — . . 


V = Via + ^ 


( 20 ) 


is used rather than the exact expression 


V = ^(Via + u)2 + v2 


( 21 ) 


The compressible velocity distribution has been computed for the cascade 
example of figure 4 by both equation (20) and equation (21) using the 
Prandtl-Glauert rule for obtaining the ccmpresslble perturbation- 
velocity components (with k = fl) . 

Uc = ui/fi Vc ^ Vi 

The resulting velocity distributions are shown in figure 5, The appre- 
ciable differences over a large part of the upper surface are a result 
of the departures of the direction of thd tangent to the upper surface 
from a direction parallel to the mean velocity. (See fig. 4.) Both 
velocity distributions have been faired over a small region near the 
nose and the tail of the blade to give a single stagnation point in 
these regions, but this fairing was over such a small region as to be 
hardly detectable on the scale of figure 5. Because of the appreciable 
difference between the two velocity distributions, a corresponding dif- 
ference in the value of the eiror expressions might be expected. ■ The 
values of the error expressions using equations (14) and (15) with 
y = 1.4 and based both on the exact equations for the resultant veloc- 
ity and the velocity components and on_the corresponding equations 
involving only the linear terms in the pertiurbation-velocity components 
are given in the following table: 


Method 

Error 

(percent) 

5m 


5x 

5r 

Linearized velocity relation 
Exact velocity relation 

-0.11 

0.02 

-15.34 
- 6.63 

-12.23 
- 3.72 

-18.91 
- 4.45 
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The Prandtl-Glauert riile, with the use of the linearized relation 
for the velocity, gives error expressions, of which the maximum for this 
example is of ahont the same order of magnitude as for the Inconpressihle 
approximation to compresslhle flow. Thus, apparently little is to he 
gained in the acctiracy of the velocity distribution by the use of the 
Prandtl-Glauert irule in this form when applied to typical compressor 
blading of moderate camber. The results slitjuld be even less accurate 
for typical turbine blading in which the camber is generally much larger. 

The reason for this poor performance is clarified by the comparison 
of the inconpressihle velocity distribution based on the same linearized 
approximation with the original inconpressihle velocity distribution in 
figure 6. The difference is of the same order of magnitude as that 
between the two interpretations of the Prandtl-Glauert rule shown in 
figure 5. The values of the error expressions for the linearized 
approximation to the inconpressihle flow, as compared with those for the 
theoretically exact incompressible flow (which represent computational 
inaccuracies) are as follows: 


Method 

Error 

(percent) 

5m 


Sx 

Sp 

Linearized velocity relation 
Exact velocity relation 

0.26 

0.00 

-9.77 

-0.09 

-11.92 
- 0.06 

-17.90 
- 0.07 


The preceding table together with figures 5 and 6 indicates that the 
errors due to linearization are quite large. In fact, these errors are 
of the same order of magnitude as the compressibility correction for 
this cascade example as is seen by comparing the linearization error 
shown in figure 6 with the conpressibillty correction shown in figure 7. 
The considerably better restilts obtained from the Prandtl-Glauert rule 
when using the exact relations than when using the linearized relations 
is thus a result of elimination of the direct error due to linearization 
for the incompressible flow. There remains, however, the error due to 
linearization in conputlng the compressibility correction. Whenever the 
direct errors due to linearization are as large as or larger than the 
compressibility correction, more accurate results may be expected if 
exact relations for computing the resultant velocities and velocity com- 
ponents from the perturbation velocities are used. 

Approximation based on linear pressure-volume relation . - As a 
third type of approximation method to be investigated, the inverse 
method of reference 3 based on the use of a linear relation between 
pressure and specific volume was used. In this method, a proportionality 
constant is determined between the dimensionless compressible velocities 
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(ratios of fluid velocity to stagnation velocity of sound) for the 
approximated compressible flow of a polytropic gas and a theoretically 
exact compressible flow of an ideal fluid with a linear relation between 
pressure and specific volume. The proportionality constant is deter- 
mined in such a way that the continuity equation is satisfied for the 
approximated flow, that is, = 0. Because the flow with the linear 
pressure-volume relation is ii^rotational, the approximated flow will 
also be Irrotational and thus satisfy the circulation condition (4), 
that is, = 0, within the numerical accuracy of the solution. The 

momentum conditions (2) and (3) will not, however, be exactly satisfied 
and the values of — 8y and 6^ will therefore serve as a measure of the 

acctiracy of the approximation. 

In order that this approximation method might reasonably be com- 
pared with the previously considered methods, it was desirable to obtain 
a cascade geometry and flow canditlons similar to the other approxima- 
tions. Because this approximation as well as the inconipressible-flow 
approximation was obtained from inverse design methods, it was impossible 
to obtain exactly the same blade geometry. Furthermore, the Inconpres- 
sible approximation and the Prandtl-Glauert approximation are for some- 
what different configurations because of the different values of Sy. 

The Mach nuinber based oh the vector mean of the upstream and downstream 
velocities was chosen the same for all the different approximtions, and 
the velocity distribution for the present approximation (fig. 8) was 
chosen to be veiy nearly the same as for the example of the Prandtl- 
Glauert approximation in which the exact velocity relation (fig. 5) was 
used, but blade shape, blade spacing, and flow angles are somewhat dif- 
ferent (figs. 4 and 9). However, the differences between the different 
approximation examples are small enough to make relative values of the 
error expression signlflcant^in comparing the relative accuracy of the 
different approximation methods. 

As in the case of the inconpressible-flow solution, it is desirable 
to have some measure of the numerical inaccuracies of a particular solu- 
tion Independent of the errors in approximating the desired flow of the 
usual polytropic gas (y = 1.4). The numerlpal Inaccuracies in obtain- 
ing the flow with a linear pressure-volume relation (y = -l) can 
readily be obtained from the general error expressions (lO) to (l3) by 
using the proper equations for dimensionless pressure and density for a 
gas with linear pressure-volume relation 


p = A - yiTqS 


1 + 


( 22 ) 


P 


(23) 
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where A is an arhltrary constant and q. is the dimensionless velocity 
(ratio of velocity to stagnation velocity of sound) for this gas. The 
constant A does not affect the values of any of the error expressions 
for closed "blade profiles and was therefore taken as zero. 


For the present cascade example, the error expressions have "been 
calculated both for the linearized pressure-volume flow (r = -1) -by 
using equations (22) and (23) and for a polytropic gas with y = 1.4 by 
using equations (l4) and (is); the results are shown in the following 
table: 


Method 


Error 


(percent) 



6m 


5x 

Sp 

Linearized pres sure -volume 
Approximated flow for y = 

flow 

1.4 

0.00 

0.00 

0.47 

-0.33 

0.46 

-0.99 

0.00 

0.00 


Because this method is theoretically exact for the linearized pressure- 
volume flow, the error expressions for this case are measures of the 
numerical inaccuracies due to rounded-off values, integration approxi- 
mations, etc. The error expressions for the other case (y = 1.4) are 
a measure of the combined error due to numerical Inaccuracies and the 
inherent error resulting from assuming the velocity Vg for this flow 
to be proportional to the velocity for the linearized pressure-volume 
flow (y = -l) . The errors for the approximated flow are apparently 
or0.y slightly greater than the numerical inaccuracies in the calculation. 
The error expressions have also been calculated for all the examples 
given in reference 3 and for all cases were below 1 percent. Thus, 
apparently this method is quite accurate for a wide variety of examples 
and is not limited to flow with small perturbations from the mean veloc- 
ity, as is the case with the Prandtl-Glauert approximation. 

For convenient corrparison, the error expressions for the examples 
of the different approximations to flow with y = 1.4 previously pre- 
sented are summarized as follows: 


Method 


Error 

(percent) 





5x 

8p 

Inconpresslble approximation 
Prandtl-Glauert, linear 

6.72 

2.05 

-14.14 

- 0.07 

velocity relation 
Prandtl-Glauert, exact 

-0.11 

-15.34 

-12.23 

-18.91 

velocity relation 
Approximation based on linear 

0.02 

- 6.63 

- 3.72 

- 4.45 

pressure -volume relation 

0.00 

- 0.33 

- 0.99 

0.00 
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DISCUSSIOT 

It should be emphasized that the error expressions measure how 
closely certain necessary conditions, representative of particularly 
significant features of the flow, are satisfied; but even if they are all 
zero, the solution may still be inaccurate. An Infinite numoer of con- 
ditions are required for a solution to be exact; namely, that it agree 
with an exact solution everywhere. Inasmuch as the error expressions 
are based cn integrals, they are essentially particular methods of 
averaging the errors of a solution. In certain pathological cases, the 
errors may cancel each other in the Integration, giving a null value for 
the error expressions. If they are used with discretion, however, the 
error expressions furnish valuable checks on the accuracy of compressible- 
flow solutions. 

Althou^ the expressions may be used in checking any type of 
approximate solution, they are particularly valuable in checking solu- 
tions Trfalch give the velocity or pressure distribution on the blade sur- 
face directly, as there are often no other convenient checks available 
in this case. For relaxation solutions or other solutions which give 
the flow throughout the field, the residual of the difference equation 
approximating the differential eqxiaticn furnishes other indications of 
the error, but does not indicate the errox involved in approximating the 
differential equation. Hie error expressions thus furnish valuable 
supplementary checks in these cases. 

For convenience of discussion, it has previously been assumed that 
the approximate flow solution is given in the form of the dimensionless 
velocity distribution, but this is unnecessary. Hie solution may be 
specified in any one of a variety of ways, such as the upstream and 
downstream Mach numbers and a pressure-coefficient distribution around 
the blade. The error expressions would then be calculated using the 
relations for density, pressure, and velocity in terms of the specified 
variables which are exact for the fluid which it is desired to approxi- 
mate. 


Lewis Flight Propulslcn Laboratory, 

National Advisory Committee for Aeronautics, 
Cleveland, Ohio, June 28, 1951. 
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APB3NDIX A 
SYMBOLS 

The following symbols axe used in. this jreport: 

A arhitraxy constant 
a velocity of sound 

B Integration path along hlade contoux 

C Integration path enclosing only fluid 

k arbitrary constant 

M Mach nimber 

p pressure 

q ratio of velocity to stagnation velocity of sound for gas with 
linear pres sure -volume relation 

S cascade spacing 

s arc length 

u pertvixbation-velocity component parallel to mean velocity 

V velocity 

V perturbation-velocity conponent normal to mean velocity 

X,Y rectangular coordinates, parallel and normal, respectively, to 
mean velocity 

AY-g range of Integration of Y on blade 

x,y rectangular coordinates normal and parallel, respectively, to 
cascade axis 

X ratio of specific heats 

error in mass flow, percent (equation (lO)) 

error in x-can 5 )onent of blade force, percent of resultant force 
(equation (12)) 
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error in y-coniponent of "blade for.ce, percent of resultant force 
(equation (ll)) 

error in circulation, percent (equation (13)) 
circle angle 

angle "between velocity vector and normal to cascade 
density 
solidity 

perturbation-velocity potential 
constant, //l - 

Su"b scripts: 

1 far upstream of cascade 

2 far downstream of cascade 

B along blade contour 

c compressible 

i incompressible 

m values with, reference to vector mean of upstream and downstream 
velocity 

mB.y maximum value 

s conrponent in s-direction 

T total or stagnation conditions ■ 

X camponent in X-dlrectlon 

X component in x-dlrection 

Y conponent in Y- direction 

y component in y-dlrection 


6p 

0 

A 

P 

a 

Q 
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APEE2IDIX B 


EROOF THAT EQUATIONS (l) TO (4) ARE SATISFIED WITHIN 
LINEAE APPROXIMATION FOR HIANDTL-GLAUEBT RULE 

It •will be sh.o’wn that equations (l) to (4), ■which, are evaluated 
along contour C of figure 1, are sa'bisfied for the Ptandtl-Glauert 
approximation within linear terms of the perturbation-velocity com- 
ponents} that is, they are satisfied to the approximation that higher- 
order terms in the con^jressible and incompressible perturbation-velocity 
components may be neglected in comparison -with the linear terms in com- 
puting pressures, densities, and velocities for the compressible and 
Incon^iressible flows. Equa'tlons which are accurate only within linear 
terms in the perturbation velocity will be indicated by the sign « . 

In establishing that these equations are satisfied within this 
linear approximation, the equations are for convenience expressed in 
terms of the coordinates X,Y paral l el and normal, respectively, to the 
mean velocity rather than in terms of the coordinates x,y normal and 
parallel, respectively, to the cascade axis. The expressions in the new 
coordinates may readily be obtained from the general equations (l) to (4) 
in the same manner as equations (5) to (8) were obtained. The resulting 
expressions are 


- (pVx)2 - (pVx)l Sy - (pVy)2 - (pVy)^ = 0 (Bl) 


P dX - (pg- Pi)Sx - [(pVxVy)2 - (pVxVi]Sy 

- [(pVy^)2 - (pV^I ° 


p dY - (pg- Pi)Sy - |(pVx^)2 ■ ^ 

- [(pVx's - ° 
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V, ds + = 0 


In order to check whether these equatlonB are satisfied within the 
linear terms in the perturbation-velocity components u,v defined as 



u = 


Vx - 


'm 


and 


V = Vy 

the equation will now be expressed in terms of these pertubatlon-velocity 
components in which only linear terms will be considered. The values of 
^m^ Pjjjj Sy, S, and chord are taken as standard or zeroth order. 

The perturbation-velocity components u and v for both the compres- 
sible and incompressible flow are assumed to be of the first order. 
Camber^ thickness, range of Y on a blade, and turning angle must then 
be of the first order because of the relation along any streamline 

„ JL 

The linearized expressions are obtained as follows: 



~1 . Ili _z!" 

r-1 

1 . Ill !s! /i + _ I::! Ji£ 
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2 „ 2 
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Pm 

1 Izi, 
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o o 




2 Bqi2 


(r-DVm^ u' 


JL. 

ir-i 


1 - 




'm 


1 
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(r-i)v^2 


JL. 

ir-i 








(B5) 
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Hence 


P-Pm « 


Vja 


Similarly 



- iv" ^ 



« P 


V 2 
m m 



p V ^ 


1 + ( 2 -\ 2 ) 



pVxVyfJpj/jji^ ^ 

m 



ID. 


'm> 


K 0 



r 1 


r- 


- ~ 


(^P dj 

fS 


< 



.Jb 


Jb 


_ 


(B6) 

(B7) 

(B8) 

(B9) 

(BIO) 

(Bll) 

iB12) 
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where AYj indicates range of integration of Y on the hlade. Since 
"both Alg and are of the first order, the integral is equal to 

zero to the first order. This merely means that the blade force is 
normal to the mean velocity within the linear approximation, a condition 
which is exactly satisfied for incompressible flow. If these values are 
substituted in equations (Bl) to (B4), the linearized form of each of 
these equations for conipresslble flow results. 



Because equations (B13) and (B15) are identical and equations (B14) and 
(B16) are Identical, except for factors involving and Vjj^, which 

are of standard order, the four independent-conditions for the general 
case (equations (Bl) to (B4)) are reduced to two Independent conditions 
for the linear approximation. 

In proving that these conditions are satisfied for solutions 
obtained by the Erandtl-Glauert rule, they are expressed as equivalent 
conditions on the incompressible flow by the Erandtl-Glauert transforma- 
tion and it is then proved that these resulting conditions are satisfied 
within the linear approximation for the exact solution of the incompres- 
sible flow. 


2251 
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From the haBic Prandtl-Glauert relations (equation (l9))^ the fol- 
lowing additional relations between the corresponding compressible and 
incompressible flows are obtained (reference 2): 


Sx 




u 


c 


= Sx,i 

%,c = 

(B17) 

k ^ 

^c = § ^i 

(B18) 


The substitution of these values in the two Independent equations (B13) 
and (B14) gives the conditions on the compressible flow expressed in 
terms of the corresponding inconpressible flow as follows: 




'U2-U1' 


'm 


I Sy . 
'i 


V2-V1 


V, 


m 




« 0 


(B19) 


Pm^m 


2 1 



U2-U1' 


'm 


%,i 


V2-V1 


'm 




0 


(B20) 


The incompressible flow on which the Prandtl-Glauert approximation is 
based must satisfy the conditions (Bl) to (B4) with p constant and p 
given by equation (I8) . As in the case of the compressible flow, the 
linear approximation to these conditions gives only two Independent con- 
ditions, which may be most easily obtained by considering these condi- 
tions as the limiting case of a compressible flow as the Mach number 
approaches zero in equations (B13) and (B14) . Thus 


PmVi 


/U2-U1 


m'm 


'm 


{ - (' 


V2-V1 


'm 




S! 0 


(B21) 


PmV^ 
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0 If). ‘^B,i - 

L^B,i ^ ^ 


U2-ui\ 


V, 


m / X 




a 0 (B22) 
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Equation (B19) differs from equation (B2l) by the factor i^/k and 
eq-uation (B20) differs from equation (B22) by the factor l/k. Thus, 
if Ci/k and l/k are not large, the conditions on the compresslble- 
flow equations (B19) and (B20) obtained from the Prandtl-Glauert rule 
are satisfied to the same order as the linearized conditions on the cor- 
responding incompressible flow. According to equations (B18), this con- 
dition is equivalent to the condition that the compre-ssible 
perturbation-velocity ccn^jonents Uq and v^ are of the same order of 
magnitude as the corresponding components for the incompressible flow. 

For the linearized form of the necessary conditions on the incom- 
pressible flow to be satisfied within a given accuracy, the incompres- 
sible perturbation-velocity^ comqjonents must be sufficiently small com- 
pared with Vjn and, consequently, blade-surface slope, blade camber, 
blade thickness, angle of attack, turning angle, etc. of the incompress- 
ible cascade must be sufficiently small. In order to guarantee that 
the corresponding compressible flow satisfy the necessary condition with 
similar accuracy, the compressible perturbation velocities u^ and v^ 
must also be sufficiently small. If, in a given Incompressible cascade, 
the values of ui and v^ are equal to the maximum allowable values t 

for the required accuracy and it is required that values of u^ and v^, 
for the corresponding compressible cascade not exceed these values, it 
is necessary that k > 1. Inasmuch- as the ratio of the compressible- to 
the incompressible values of blade camber, blade thickness, angle of 
attack, turning angle, etc. are equal to 0/k, the ratio of these quan- 
tities must not be greater than ^l-Mja^ and, consequently, blade thick- 
ness, turning angle, etc. for the compressible flow must be extremely 
small if % 1 order for the Prandtl-Glauert approximation 

to be reasonably accurate. 
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Figure 1. - Path of integration through cascade of blades used in 

general analysis. 
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Figure 2. - Velocity dletrltutlon for iucoinpresslble flow and inconrpreBsl'ble approximation. 
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Figure 6. - Veloeity distributions for exact Incompressible flow and for linearized Incompressible 
flow, showing error resulting from linearization. 
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Figure 8. - Velocity distri'bution for approximation with dimensionless velocity proportional to 

that of llnearized-pressure-volimie flow. 
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Figure 9. - Cascade geometry for approx 
prpportiorial to that of llnearlzed-pr 
solidity, 0.9856. 
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